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'^ ! Abstract 

O ' Using a bit string model, we show that asexual reproduction for 

diploids is more efficient than for haploids: it improves genetic ma- 
terial producing new individuals with less deleterious mutations. We 
also see that in a system where competition is present, diploids dom- 

Q"^ ■ inate, even though we consider some dominant loci. 

^. 
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1 Introduction 



In a system where population is kept constant, competiton just allows the 

. best individuals to survive, i.e. individuals which have less deleterious muta- 

'"O ' 

jH ■ tions and are better adapted to the current enviroment have a higher proba- 

cy . bility to survive. Haploid gametes reproduce by mitosis, mutations being the 

only source of variability, thus they cannot improve too much when repro- 
ducing. On the other hand, diploids could reproduce with crossings between 
rS I two equivalent copies of genetic material, by meiotic parthenogenesis. This 

c^ • is supposed to be a more efficient way to improve the ability to adapt the 

population to the environment. The purpose of this work is to verify this 
hypothesis within a simple bit string model [Q of population dynamics. 

The hypothesis of Jan, Stauffer and Moseley ( JSM) proposes that hap- 
loid populations convert ffist into asexual diploid and then diploids begin 
to reproduce in a sexual way. This complete process has been studied in 



references p, Q, although they used a completely different model, with con- 
versions (haploid — > diploid and asexual —>■ sexual) happening only for indi- 
viduals with many deleterious mutations. In this case, sexual diploids finally 
dominate the system with a small fraction of asexual diploids and haploids. 
In this work we show that in a system under competition, with an asexual 
population completely haploid at the beginning, and where with a small 
probabihty two haploids can fuse to create a diploid, evolution leads to a 
population dominated by diploids at the end. We also show that even the 
introduction of some dominant loci into the genetic material of diploids does 
not change too much this final result. 

2 Models for the haploid-diploid cycle 

We represent each haploid by a bit string of 1024 "0"'s and "l"'s, with a 
probability to survive given by: pi = x^'~^^, where Ni is the number of "l"'s 
in the bit string representing individual i. Every time step each haploid 
reproduces by duplicating its bit string and some Mh random mutations are 
performed in this new individual, say Mh = 1, by flipping Mh randomly 
chosen bits from to 1 or vice-versa. 

After this breeding process applied to all individuals, in order to keep 
the population almost constant, fluctuating around some number Po, say 
Po = 10000, we determine the value of x by solving equation: 

oo 

y: i/(iv)x^+i = p, (1) 

N=0 

where H{N) is the number of haploids with A^ "l"'s into their bit string. 
Already knowing the value of x, we test each individual i of the population, 
to keep it alive according to the probability x^*"*"^, killing it otherwise. All 
this represent one time step, where all individuals reproduce first and then 
each one is tested to survive (including newborns). 

For a population of diploids, we represent each individual by two parallel 
bit strings of 1024 "0" 's and "1" 's. At each time step, every diploid reproduce 
by meiotic parthenogenesis. Both strings are cut in the same random position 
and a new bit string is constructed by joinning one randomly chosen piece 
with the complementary piece of the other string. Then, we copy this new bit 
string in order to have two bit strings representing this new diploid. Finally, 
Md mutations, say Mo = 2, are randomly performed. 

The probability to survive in this case is also given by: pi = x^'^^, but 
now Ni represents the number of coincident positions with "l"'s in both bit 
strings (homozygous loci). We can also consider some dominant positions 



along the bit strings: in this case, just one 
enough to contribute to A^j. 



bit at this position (locus) is 



2.1 Populations without exchange 

First we study how two populations (one haploid and the other diploid) de- 
velop living together under the same environment, but with no possibility 
of haploids becoming diploids or vice-versa. We begin a population of Po/2 
haploids and Po/2 diploids. At each time step, both kinds of individuals 
reproduce following the rules already quoted, using the same x, and we per- 
form the sum of equation ([|) over all haploids and diploids in order to keep 
the whole population constant. 
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Figure 1: Haploids versus Diploids 

Figure 1 shows the number of haploids and diploids evolving with time. 
Most of the time diploids dominate the whole system, and eventually no 
haploid survives. The same behaviour was observed for different runs. 

This situation does not change too much if we include some dominant 
positions for diploids, that is, less than 50% of dominant positions. Only for 
an unrealistic situation where more than 50% of positions are dominant, we 
can see some instances where the haploids dominate the system. 



2.2 Exchange (diploids with only recessive positions) 



2.2.1 No return 

Now we begin with a population totally haploid, but each individual has a 
small probability Phd, say Phd = 0.01, to randomly choose another individ- 
ual to fuse and form a diploid, instead of normal reproduction. The fused 
diploid has the two bit strings of those haploids. Otherwise, this haploid 
reproduces in the normal way. During the evolution, diploids also reproduce 
according to the same rules explained before. 

In this case we can see that the population becomes totally diploid, after 
few steps. Using different values of Phd we just change the number of steps 
needed to converge to an only-diploid population. 

2.2.2 With return 

Now, we add some probability Pdh, that a diploid transform itself into two 
haploids, each one with one of the bit strings from the diploid, with some 
Mh mutations. 

In Nature this backward transformation is supposed to be harder to hap- 
pen than the first case (haploid -|- haploid — ^ diploid). The diploid organism 
needs to divide its genetic material in just two equal parts, otherwise at least 
the smallest part which misses genetic material will die before reproduction. 
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Figure 2: Fraction of diploids in population as a function of Pdh 



For this case we observed that diploids still dominate the system but there 
is a small fraction of haploids at the end, which are products of the continuous 
transformation of some diploids into haploids. In figure 2 we see the fraction 



of diploids in the population as a function of the back probability Pdh- These 
results are obtained as an average over many steps, after equilibrium has been 
reached. Althougth some extreme values of Pdh are taken (as Pdh = Phd 
corresponding to the rightmost point), we can see that diploid individuals 
always dominate the population. 

2.3 Exchange (diploids with recessive and dominant 
positions) 

2.3.1 No return 

Now, we repeat the same experiment explained before, but consider some 
percentage D of dominant loci for the diploids. As before, some haploids 
become diploids at each step, with no diploid becoming haploid. As a re- 
sult, again, after some steps the whole population becomes diploid with no 
haploids in the final equilibrium configuration. This result is independent of 
the value of D: even an extreme value oi D {D = 100%) just delay the final 
situation by some steps. 

2.3.2 With return 

If we let diploids to return back to haploids with some probability Pdh and 
we also consider some dominant positions for diploids, the results obtained 
before do not change qualitatively, as we can see in figure 3. The fraction of 
diploids at the equilibrium situation is plotted for different values of D, and 
for fixed values of Phd = 0.01 and Pdh = 0.005, taken as an example. 

3 Discussion 

In order to study the most efficient kind of behaviour (haploid versus diploid), 
we considered systems where competiton allows to survive only a pre-determi- 
ned number of individuals. Of course, these selected individuals will be the 
best adapted to the environment, presenting a small number of deleterious 
mutations. 

Diploid reproduction is more efficient to improve new generations than 
haploid reproduction. The final population is dominated by diploids, and 
haploids are present only if we let some diploids to become haploids. 

The presence of some dominant loci in diploids does not change the re- 
sults. In the limiting case Pdh = (no diploid returns back to haploid 
condition), dominance just modifies a little bit the equilibration time, but 
eventually we have no haploids anyway. 
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Figure 3: Fraction of diploids in population as a function of D 



